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Abstract. The Beta Beam, the concept of generating a pure and intense (anti) neutrino beam by letting accelerated radioactive 
ions beta decay in a storage ring called the Decay Ring (DR), is the basis of one of the proposed next generation neutrino 
oscillation facilities, necessary for a complete study of the neutrino oscillation parameter space. Sensitivities of the unknown 
neutrino oscillation parameters depend on the DR's ion intensity and of its duty factor (the filled ratio of the ring). Different 
methods, including analytical calculations and multiparticle tracking simulations, were used to estimate the DR's potential to 
contain enough ions in as small a part of the ring as needed for the sensitivities. Studies of transverse blow up of the beams 
due to resonance wake fields show that a very challenging upper limit of the transverse broadband impedance is required to 
avoid instabilities and beam loss. 
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INTRODUCTION 

The discovery of neutrino oscillations [1] has confirmed 
that neutrinos are massive and that their flavor (v^, 
v^, v^) and mass eigenstates (Vi, V2, V3) are mixed. 
Neutrino physics is now in an era of precision mea- 
surements of the parameters that govern these oscilla- 
tions: two Am?- = m^. parameters (|Am32|, |Am2i |), 
three mixing angles (023, ^12, ^13) and a CP violating 
phase {5cp). The most precise determinations to date 
are Am^i = (7.586+g-2i2) . ^y2^ ^ (34.1+}^ 

[2], |Am^2l = (2.43 ±0.13) • IQ-^ eV^ [3] and ^23 = 
(45 ± 3.4)° [4]. This leaves three unknowns: ^13 (<12.2° 



Assumed Fluxes: I . I e 1 8 v/year from '^Ne and 2.9e 1 8 anti-v/year from ^He 
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FIGURE 1. The ^13 (left) and 3cp (right) sensitivities for 
different suppression factors (SF) [5]. A SF of 0.01% (green) 
gives almost the same result as the case of no atmospheric back- 
ground (orange). The 0.1% SF (blue) gives similar sensitivity 
and with 1% SF (red) the sensitivity has decreased slightly. 

[6]), 5cp and sign(Am32). The discovery of ^13 > and 
^CP 7^ 0° and 7^ 180° would mean existence of CP vi- 
olation in the leptonic sector. Near future experiments 
will perform precision measurements of the knowns and 
continue to probe the unknowns (e.g. T2K started taking 
data 2009 [7]). However, a complete study of all neu- 
trino parameters requires a better characterized neutrino 



beam with higher flux then ever available before. One 
of the three present options for a next generation neu- 
trino oscillation facility [8] is the Beta Beam concept [9] 
wherein it is proposed to store high energy (7 = 100) ra- 
dioactive ions in a horse-track shaped storage ring, called 
"Decay Ring" (DR), with a straight section pointing to 
a neutrino detector. Ions that beta-decay in the straight 
section emit electron (anti) neutrinos in a pure (v^) 
beam with opening angle \/y. The aimed annual (anti) 
neutrino fluxes of (l.lel8) 2.9el8 [10] from (jS") j3+ 
decaying (^He) ^^Ne ions give ^13 and 5cp sensitivities 
shown in fig. 1 [5] for different suppression factors (SF) 
of the detector (which coincides with the duty factor of 
the DR). Fig. 1 shows that with the nominal fluxes the 
beam can only fill less than 1% of the DR. The challenges 
of how to produce enough ions, how to accelerate them 
through a Beta Beam complex and how to achieve 0.58% 
SF are described elsewhere [10]. In the studied scenario 
the Beta Beam complex is based at CERN, making use 
of already existing PS and SPS machines and the DR has 
the same circumference as SPS, C = 691 1 .6 m. Assuming 
20 bunches of 6911.6mx0.58%/20 = 2 m each, the num- 
ber of i^Ne (^He) per bunch have to be 3 • 10^^ (4 • 10^^) 
to reach the nominal (anti) neutrino fluxes. In this report 
we focus on the question whether the required amount of 
ions can populate such short bunches without too large a 
risk of beam instabilities. This collective effect study is 
based on the previous design studies [10] and all param- 
eters used are listed in table 1 and 2. 

COLLECTIVE EFFECT STUDIES 

High intensity ion beams, foreseen for the Beta Beam 
project, could suffer "Collective Effects". These are 
caused by electromagnetic interactions between particles 
in the beam, either with each other directly or through the 



TABLE 1. Input parameters (some from [10]) above 
the line and calculated parameters below the line. 



Parameters 


DR i^Ne 


DR ^He 


Z 


10 


2 


A 


18 


6 


Vrf [MV] 


11.96 


20.00 


Brest [MeV] 


16767.10 


5605.54 


Nb 


3.1e+12 


4.0e+12 


h/2 [s] 


1.67 


0.81 




3.60 


6.00 


ro [m] = rpZ^/A 
Etot [G^^y-Erest 

hm=ZeNBlTrev 

e]"" [eVs] = jP^Etot^Smax 


8.53e-18 
1676.71 
755.80 
0.22 
43.27 


1.02e-18 
560.55 
195.04 
0.06 
14.46 



TABLE 2. Input parameters (mostly from the previous Beta 
Beam Decay Ring design report [10]) above the first line. 
Assumed transversal impedance parameters between the lines. 
Calculated parameters below the last line. These parameters are 
the same for the two isotopes, ^^Ne and ^He. 



Parameters 


Description 


Value 


h 


Harmonic Number 


924 


C[m] 


Circumference 


6911.6 


P [m] 


Magnetic Radius 


155.6 


Ytr 


Gamma at Transition 


27.00 


7 


Relativistic Gamma 


100.0 


^max 


Max. Mom. Spread 


2.5e-3 


Lb [m] 


Full Bunch Length 


1.970 


Qy 


Vertical Tune 


12.16 


{P)y [m] 


Av. Ver ptron Tunc. 


173.64 


by [cm] 


Ver Beam Pipe Size 


16.0 


fix 


Quality Factor 


1.00 


(Or_L [GHz] 


Ang. Resonance Freq. 


6.28 




Shunt Impedance 


2.00 



Tl = rtr -7 

Trev [M s]=C/(/3c) 


Phase Slip Factor 
Revolution Time 


1.27e-3 
23.0558 


(Orev [MHz] = ^ 


Ang. Revolution Freq. 


0.27 


^ / /zZeV 7] COS 05 


Synchrotron Tune 


3.63e-3 


(Oc [GHz] = 1^ 


Cut-Off Ang. Frequency 


1.87 



surrounding environment. Collective effects could limit 
the final performance of the accelerators. The studies of 
instabilities of all ion beams and all machines in the Beta 
Beam complex is therefore a crucial part of the project. 
Here we have focused on ^^Ne and ^He in the DR. 

A particle traveling inside an accelerator leaves elec- 
tromagnetic fields lagging behind. Trailing particles feel 
a force due to the net field caused by all preceding parti- 
cles. The line integral of this force over a certain length 
(which could be a part of the beam chamber) gives what 
the particles see as "wake fields". If the wake fields last 
for the duration of the bunch (^100 ps) particles in the 



"tail" of the bunch can interact with the wake fields 
caused by the particles in the "head" of the bunch and 
cause single bunch instabilities. 

The action of the wake fields are described by the 
wake potential, W{t), in the time domain and by the 
impedance, Z{co) = ^ [W{t)], in the frequency domain. 
This report studies impedances caused by wake fields 
trapped in cavities of the vacuum chamber, so called 
resonance impedances, Z^^^{(0). If the quality factor is 
Q = Ry/C/L and the resonance frequency is 0)^ = 1/ V^LC 
the resonance impedance can be modelled as an RLC 
circuit [1 1] in the transverse plane as 
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where R± is the transverse shunt impedance, assumed to 
have a value close to RHIC; 2 MQ/m [12]. So far we have 
only studied short lived resonance wake fields, i.e. broad- 
band (2=1) impedances in the transverse plane. There 
are many different types of collective effects that could 
lead to beam instability but this study is constrained to 
transverse broadband resonance impedances. Three dif- 
ferent methods have been used to achieve the maximum 
number of allowed ions per bunch, N^^. 

One approach is to use the peak current values of the 
bunch current and momentum spread as input to a coast- 
ing beam formula. This gives an expression for the in- 
tensity limit that we will call the coasting beam equa- 
tion [13] and which for zero chromaticity (as assumed 
all through this report) is 
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Here c is the speed of light in vacuum and all other 
parameters are given in table 1 and 2. 

IMOSES [14] solves an integral equation in the fre- 
quency domain to give the rise time, T, of the instabil- 
ities for different head-tail modes as a function of the 
bunch intensity. The limit, /^^, is given by the most cru- 
cial head-tail mode after defining the maximum allowed 
growth rate, ( 1 / t)^^ . To reach the ion equivalent intensity 
threshold we divide by a factor Z; = /f /Z. The maxi- 
mum allowed number of ions per bunch is then given by 
the conversion A/g^ = TrevP^ jZe. 

The third method uses the multiparticle tracking code 
HEADTAIL[15] where a bunch of macro-particles is 
sliced longitudinally and the impedance is assumed to 
be localized at a few positions around the ring. At each 
impedance location, each slice leaves a wake field be- 
hind and gets a kick by the field generated by the pre- 
ceding slices. The bunch is then transferred to the next 
impedance location via a transport matrix. For the Beta 
Beam studies the possibility of bunches with ^^Ne and 



^He was added to the code. An exponential least squared 
fit to the envelope of the vertical oscillation of the mean 
bunch position gives the growth rate of the instability. 
Same as for MOSES the bunch intensity limit, N^^, is 
reached when the rise time is shorter than allowed, i.e. 

It could be argued that instabilities with the longest 
rise times should define A/^^, i.e. 0. However 

in this report we will take an optimistic approach, assume 
that slow instabilities can be damped with sextu- and oc- 
tupoles and define = 400 1/^ for both MOSES 

and HEADTAIL. 

With the three methods, mentioned above, we studied 
the effect on the bunch intensity limit, A/^^, by changing 
slightly the longitudinal emittance, £/, (fig. 2 (a) and 

(b) ) and assuming = 2 MQ,/m. Fig. 2 (a) shows that 
according to MOSES (HEADTAIL) increasing £/ with 
about 5 (10) eVs from the working point (indicated by 
grey arrow) the desired number ^He per bunch, 4 • 10^^, 
would be acceptable. This would however also mean an 
undesired increase in SF and momentum spread (also 
indicated in fig. 2 (a)). As can be seen in fig. 2 (b) the 
bunch intensity limit for ^^Ne, 3 • 10^^, is far out of reach 
when = 2 M^l/m is assumed. 

Since impedance could improve in modem machines 
compared to old accelerators a scan over the shunt 
impedance was performed to see the impact on N^^ (fig. 2 

(c) and (d)). Fig. 2 (c) shows that for a shunt impedance 
at the level of SPS, R'f = 20 M^l/m, maximum number 
^He allowed per bunch, according to HEADTAIL and 
MOSES, is not more than 300- 10^ For A^^ = 4 • 10^^ ^He 
per bunch R± ^ 2 MQ./m (similar to RHIC) is needed. 
However, as can be seen in the log-log scale plot of fig. 2 

(d) R^ ^ 0.2 MQ/m is needed in the DR to allow 3 • 10^^ 
^^Ne per bunch. 

Scans over resonance frequency, fr = (Or/ {2%), and 
chromaticity, t,y = {dQy/Qy)/{dp/ p), were also per- 
formed without any significant relaxation in bunch in- 
tensity limit, N^^, within realistic ranges of the scan pa- 
rameters. 

CONCLUSIONS 

There will be large challenges due to requirements 
of seemingly insurmountable low transverse broadband 
impedance of the Beta Beam Decay Ring. This study, 
based on parameters mostly from [10] (table 1 and 2), 
suggests a reoptimization of the Beta Beam Decay Ring 
design. 
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FIGURE 2. N^^ as a function of (a, b) the longitudinal emit- 
tance and (c, d) the transversal shunt impedance according to 
the C.B. eq. (2), MOSES [14] and HEADTAIL [15]. The right 
(left) column shows the case for the (anti) neutrino emitter, ^He 
(^^Ne), in the Decay Ring. With the log-log scale in (d) we see 
that Rj^ ^ 0.2 MQ./m is needed to allow 3 • 10^^ ^^Ne per bunch. 
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